The late-time tails in the Reissner-Nordstrom spacetime revisited 
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We propose that the late-time tail problem in the Reissner-Nordstrom (RN) spacetime is dual to 
a tail problem in the Schwarzschild spacetime with a different initial data set: at a fixed observation 
point the asymptotic decay rate of the fields are equal. This duality is used to find the decay rate 
for tails in RN. This decay rate is exactly as in Schwarzschild, including the case of the extremely- 
charged RN spacetime (ERN). The only case where any deviation from the Schwarzschild decay 
rate is found is the case of the tails along the event horizon of an ERN spacetime, where the decay 
rate is the same as at future null infinity. As observed at a fixed location, the decay rate in ERN is 
the same as in Schwarzschild. We verify these expectations with numerical simulations. 

PACS numbers: 04.70Bw, 04.25. Nx, 04.30.Nk 



I. INTRODUCTION AND SUMMARY 

Late-time tails in black hole spacetimes have been 
studied extensively since Price's seminal analysis 
Complete understanding is available for the 
Schwarzschild spacetime — where it was also found in 
fully nonlinear numerical simulations @, HJ — and much 
progress has been made recently also in understand- 
ing the tail problem for Kerr black holes 0, and 
even for spacetimes without any assumption of symme- 
try (but with a globally weak field) [g|. The tail of 
spherically-symmetric, electrically-charged, static black 
holes — namely Reissner-Nordstrom (RN) black holes 
— was studied first by Bicak for the case of scalar per- 
turbations (for both the non-extreme and extreme cases) 
in 0] and for the case of coupled electromagnetic and 
gravitational (linearized) perturbation in Q (see also H) 
and revisited by Gundlach, Price, and Pullin (GPP) [ig |. 
Price's law has been proved rigorously in the Schwarza- 
schild case and a self-gravitating scalar field by Dafermos 
and Rodnianski (Tlj . 

The late-time tails in black hole spacetimes are under- 
stood as a result of the asymptotic form of the effective 
potential at great distances from the black hole, and the 
exponential drop off of the potential close to the event 
horizon. It is the scattering of the waves off the effective 
potential at great distances that is responsible for the cre- 
ation of the tails. (This picture applies also for the case 
of a Kerr spacetime, but this simple behavior is masked 
by an intricate mode coupling effect. See Q.) Indeed, 
it was implied first by Bicak in Ref. Q and then argued 
explicitly in greater detail by GPP in Ref. [lOJ that be- 
cause in RN the effective potential has the same form at 
great distances as in Schwarzschild, the decay rate of the 
tails in RN must be the same as in Schwarzschild. It was 
further demonstrated numerically in Ref. [l(| that this is 
indeed the case. (This argument is further reinforced by 
the analysis in Ref. [II].) 

The GPP argument suggests that all spacetimes that 
share the same asymptotic form of the effective potential, 



have the same decay rate for the late-time tails. The 
possible exception of the late-time tails in the extreme 
RN (ERN) spacetime is therefore intriguing. Specifically, 
Bicak argued in Ref. [7| that the exponent of the tails (as 
observed in a fixed observation point) in the ERN space- 
time for initial data with an initial static moment I was 
I + 2, whereas in Schwarzschild and (non-extreme) RN 
it is 21 + 2. The asymptotic form of the effective po- 
tential is the same in those three spacetimes. (They dif- 
fer only to 0(r~ 3 ), r* being the tortoise coordinate. To 
0(ln(r*/M) r~ 3 ) the three effective potentials are iden- 
tical). Applied naively, the GPP argument suggests that 
the same exponent for the tails should be observed for 
all three cases. 

However, it turns out that the effective potential in 
ERN is very different from the effective potentials in 
Schwarzschild or (non-extreme) RN near the black hole: 
instead of dropping off exponentially (in r* , for large and 
negative values of r„) towards the event horizon, the ef- 
fective potential in ERN is effectively centrifugal asymp- 
totically close to the event horizon. Moreover, the de- 
viations from centrifugality have the same leading-order 
form as the deviations from centrifugality at great dis- 
tances from the event horizon (large and positive r*). 
That is, the effective potential is asymptotically sym- 
metric 0, such that the tails at r = const are expected 
to have contributions both from scatterings at great dis- 
tances, and from scatterings very close to the event hori- 
zon. (In that sense, the event horizon is equivalent to fu- 
ture null infinity. In fact, the close analogy of the event 
horizon to future null infinity is even deeper: the tail 
along the event horizon of ERN turns out to be the same 
as along future null infinity (see below).) 

Can the contributions to the tails coming from the 
close vicinity of the event horizon in ERN overwhelm the 
contributions to the tails coming from great distances, or 
interfere with them to create tails with the indices pre- 
dicted in In this paper we study this issue, and show 
that with adaptations, the GPP argument is applicable 
also to the ERN spacetime. We then present also nu- 
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merical simulations of tails in ERN. These simulations 
indicate that the fall-off of these tails in the ERN case 
is the same as in the non-extreme case, i.e., quicker than 
that predicted in Ref. 0]. 

The organization of this paper is as follows: In Sec- 
tion [IT] we describe the scattering problem, and discuss 
the asymptotic symmetry of the effective potential. This 
symmetry was found in 7], and is presented here for 
completeness. In Section IIIII we argue that in both 
Schwarzschild and ERN the tails along r = const re- 
sult from the secondary waves (first scattering). Tertiary 
waves (second scattering) or higher-order waves that are 
present at late times originate only from scatterings at 
small values of \r*\; scattering at great distances leads 
only to secondary waves. As it is the scattering at great 
distances which leads to the tails, we argue that it is the 
secondary waves which are responsible to the tails. In 
Section ITVl we discuss the tails in ERN when the initial 
data have compact support. We argue that the tail prob- 
lem along r = const in ERN is dual to a tail problem in 
Schwarzschild (for a different, but related, initial value 
problem); as the latter is well understood, we may pre- 
dict the power-law indices for ERN tails. The meaning 
of duality here is as follows: while the value of the field 
at a particular event is different in the two spacetimcs 
(because the intermediate effective potentials of the two 
dual spacetimes are not identical), the scattering dynam- 
ics from regions of spacetime that contribute to the for- 
mation of the tail is similar, so that the same power-law 
indices for the tail are expected. We emphasize that the 
duality argument pertains only to the tail at a fixed ob- 
servation point (along r — const) . Then, in Section |V] we 
consider the tails in ERN with an initially static moment 
for the initial data. We again show that this problem is 
dual to another tail problem in Schwarzschild, and use 
this duality to find the decay rate of the tails in ERN. 
Finally, in Section IVII we present numerical simulations 
of tails in ERN which are in full agreement with our ex- 
pectations based on the arguments brought is Sections 
IIVI and El We describe the numerical code and the con- 
vergence tests done in Appendix [Al 



II. THE ERN EFFECTIVE POTENTIAL 



The ERN metric is given by 



ds 2 = 
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(1) 

where r is the regular radial Schwarzschild coordinate 
defined so that spheres of radius r have surface area of 
Airr 2 , and dfl 2 is the line-element of the unit 2-sphere. 
In terms of the Regge- Wheeler 'tortoise' coordinate r*, 
defined by 



dr* _ ( M 
dr \ r 



(2) 



the wave scattering problem is governed by the wave 
equation 

where the ERN effective potential, for a multipole £, is 
given by 

(3) 
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For a scalar field cf>, the dimensionally-reduced field ip is 
defined by 4> — rip. For field spins other than s = 0, the 
effective field tfi is interpreted appropriately [1, Q . 

To find the asymptotic effective potential, integrate 
Eq. © to find 
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and choosing the integration constant, the relation be- 
tween the radial coordinates is 



R = p + 21np 

P 



(6) 



In order to find V(r*) as r* — > ±oo, consider two cases: 
Case a) r > M (p, R > 1), for which R ~ p + 21np 
(neglecting the term — 1/p in Eq. ([6])), and Case b) < 
r - M < M (0 < p < 1, R < -1, \R\ > I), for which 
R ~ 2 lnp — p~ x (neglecting the term p in Eq. ©). 



A. Case a: r > M 

In this case 

i?~p + 21np, , 
which can be solved as 



p~2w(±e R '^ 



(7) 



where W(x) is the Lambert VF-function 1 [13| . defined by 
the inverse function of W(x) e w ^ = x (i.e., the (real) 



1 For the Schwarzschild spacetime in the usual Schwarzschild co- 
ordinates, the 'tortoise' coordinate is defined by dr* = dr/(l — 
2M/r). Denning p := r/(2M) - 1 and R := r»/(2M) - 1, and 
choosing the integration constant, the two radial coordinates are 
related by R = p + In p. In terms of the Lambert W— function we 
may write p = W[exp(-R)], or 



r(r*) = 2M + 2M W (e - 2m 
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function W(x) which solves this equation). As by defini- 
tion p > 0, the principal branch of the VF-function (de- 
noted Wo(x)) is chosen, and one does not have branch 
ambiguity. For x ^> 1, W{x) ~ In a: — In In a; + • • •, so 
that 



1 



J? — 2 In 



,R/2 



R 



2 In 



1 



,R/2 



r(r*) = r* - 2Af In H . 

Substituting in Eq. ©, one finds (r* 3> Af): 



2M 



(8) 



(9) 



V(r(n)) = + 4M^-^ In ^ + 0(r,~ 3 ) . 



(10) 



Combining the results for Cases a) and b), one finds 
(|r*| > Af) 



r? + IrJ 3 2Af + l * J 



V(r(r.)) = 

' * i'*r 

(14) 

The last expressions shows explicitly that the effective 
potential is asymptotically symmetrical under r* <-» — r*. 
Notably, the asymptotic effective potential for ERN coin- 
cides with the Schwarzschild asymptotic effective poten- 
tial (for r* 3> M). The same effective potential is also 
given in (7| . We add that in terms of the Schwarzschild 
radial coordinate r, for r>M, 

is the same in Schwarzschild as in ERN. The ERN and 
Schwarzschild effective potentials, for t = 1, are shown 
in Fig. Q] 



B. Case b: 0< r - M < M 



In this case 



R h 2 hip, 

P 



which can be solved as 



p rv, e '^'Z 2 exp 



Using the definition of the Lambert IF-function, 
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Substituting in Eq. one finds (r* < -Af): 



(11) 



(12) 



^))^-4M^lnM + 0(r ,3, 



(13) 
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FIG. 1: The effective potential for I = 1 in ERN (solid curve) 
and Schwarzschild (dashed curve), as a function of r». 

The key step in finding the asymptotic effective poten- 
tial was finding r(r*). The latter may be found without 
invoking the Lambert VF-function by an iterative solu- 
tion, in which one substitutes r*(r) into the equation 
r = ?** + g(r) . After each iteration, the right hand side 
will be a combination of terms in and in r, but the lat- 
ter terms become smaller with each iteration, and may 
be neglected. 



III. THE LATE TIME TAILS RESULT FROM 
SECONDARY WAVES 

Consider an initial perturbation field of compact sup- 
port. (We assume compact support here without loss 
of generality in order to separate at late times between 
primary waves and scattered waves, as we are interested 
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here only in the latter.) The part of the radiative field at 
later times which propagates along the geometrical optics 
rays is the primary waves. These waves scatter off the ef- 
fective potential. The waves which scatter just once are 
the secondary waves, these which scatter twice are the 
tertiary waves, and so on. (See Fig. [2]) 




FIG. 2: A schematic diagram of scattered waves. The ini- 
tial pulse in this Figure is initially outgoing (and is there- 
fore shown on the incoming leg of the characteristic hyper- 
surface), with the primary waves (p) arriving to scri+ in the 
Figure. The waves which cross the worldline of an observer 
at r — const at late times are scattered waves: the secondary 
waves (s) undergo one scattering event, and quaternary waves 
(q) undergo three scattering events. This Figure only dis- 
playes the waves approaching the observer from the right, 
hence no tertiary waves are shown. The abbreviation "EH" 
stands for the event horizon. 

In this Section we show that in both the Schwarzschild 
and the ERN cases the late-time tails are generated by 
the secondary waves, and higher-order waves are irrel- 
evant for their generation. The tails in Schwarzschild 
are generated by the scattering of waves off the effective 
potential at great distances, i.e., it is the part of the ef- 
fective potential at r* 3> M which is responsible for the 
generation of the tails. In order to study which order of 
the waves (i.e., secondary or higher) is relevant for the 
generation of the tails, we write the effective potential in 
the form 



K 

where the curvature potential 

V{U)=A ^ >- \n— + 0(rl). 



(15) 



(16) 



Here, A equals unity for the Schwarzschild effective po- 
tential. (We parametrize here the curvature potential 
with A for the analysis of the order of the waves below.) 

Secondary waves are those which scatter just once. 
Therefore, secondary waves are expected to be propor- 
tional to A. Therefore, if we study the tails as we vary 
the value of A in Eq. (fT6|) , we can determine whether the 
tails are secondary waves (linear dependence on A), or 
include higher-order contributions (deviations from the 
linear dependence on A). Primary waves are indepen- 
dent of A. In practice, we write the toy potential as 



(i-^ 



2M 



4Ml(l+l) 



r < r 
r > r 



(17) 



that is, at small distances (r < ro) the potential is 
the exact Schwarzschild potential, and at great distances 
(r > ro) the potential is approximately Schwarzschild for 
A = 1, and non-Schwarzschild for A ^ 1. (Toy poten- 
tials have been used before to study tail phenomena (in 
Schwarzschild [J and in RN and ERN [7]), although not 
the same toy potential as in here.) 

We next use a numerical code to study the dependence 
of the tails on A and ro- The code we use, and conver- 
gence tests, are described in the Appendix. It is a Cauchy 
code with each cell computed by a characteristic "di- 
amond." The code is globally second-order convergent. 
We vary ro (in practice from lOOAf to 200M), and for 
each value of ro, we vary A (in practice from 0.25 to 4). 
For each simulation (characterized by specific values of 
ro and A) we verified that indeed tails are generate d ( in 
all cases with the same tail exponent; see also Ref. [14j| ) . 
and then we recorded the value of the field in the tail 
regime for a fixed value of time. In practice, we recorded 
the fields at t = 3000M at r* = 0, but the results do 
not depend on the choice of the evaluation time or eval- 
uation point. Figure [3] shows the magnitude of the tail 
at a fixed time as a function of A for different values of 
ro. It is clear from Fig. [3] that the dependence of the 
tail on A is nonlinear for any finite value of ro- However, 
by increasing ro we find that the deviation from a linear 
dependence decreases. To find the magnitude of the field 
in the limit as ro — ► oo, we apply Richardson's deferred 
approach to the limit ( "Richardson's extrapolation" ) for 
each value of A, and then fit the extrapolated data points 
to a straight line. Fig. suggests that the tail field would 
depend linearly on A in the limit as ro — > oo. In fact, 
in Fig. [3J we look at successive simulations, in which we 
study the dependence of the field on A for scatterings 
occuring at r > ro. 

Figure [4] schematically shows the spacetime diagram 
with three curves corresponding to three different values 
of ro. For the smallest value of ro, the tails, as recorded 
at the observation point, have contributions from scatter- 
ings at all values of r > ro, i.e., from regions III, II, and I 
in Fig. 2J For the intermediate value of ro the tails have 
contributions from regions II and I, and for the largest 
value of ro the contributions come only from region I [15j | . 



FIG. 3: The magnitude of the tail in Schwarzschild at a fixed 
value of time as a function of A for different values of ro- 
Same initial data and evaluation point were used in all cases. 
*: r = 100M; □: r = 125M; X: r = 150M; o: r = 
175M; A: r = 200M. The circles (o) are the Richardson 
extrapolations of the data at finite values of ro to ro — > oo, 
and the solid line is a best fit line of the extrapolated data 
(circles). The squared correlation coefficient for the solid line 
is R 2 = 0.9995. The dotted curves are 3<r confidence curves 
of the extrapolated data. 



Figure[3]suggests that when only scatterings at asymp- 
totically large values of ro are considered, the dependence 
of the tail field on A would be linear, which implies that 
only secondary waves would be present. The tails are 
known to be the outcome of scatterings off the effective 
potential at asymptotically great distances, such that our 
results imply that indeed the tails in Schwarzschild at 
asymptotically late times are caused by waves that are 
scattered just once, i.e., secondary waves. 

Figure [5] shows for the ERN spacetime the same infor- 
mation as Fig. [3] does for Schwarzschild. We infer that 
in ERN too, the late-time tails are caused by secondary 
waves. In the ERN case the toy potential was taken as 



(1-4^ 



2M 



r <r 



^+A*m&Hr./(2M)] r>r 



(18) 



For both Fig. [3] and Fig. [5] we used momentarily sta- 
tionary initial data with compact support, but similar 
results were obtained also for other choices. In practice, 
we used initial data which are non-zero only between r*i 
and r*f, where the field has the form 



[(r* - r*i)(r* - r*f)] f 



[(' 



where r* c = (r*i + r*f)/2, and r*i = 10M, r*f = 30M. 

In fact, as for the Schwarzschild curvature potential 
A = 1, it is the linearity of the field with A in the neigh- 



FIG. 4: The spacetime diagram with three curves correspond- 
ing to three different values of ro, dividing spacetime into 
three regions I, II, and III. See the text for more information. 




FIG. 5: Same as Fig. [3] for the ERN spacetime. The squared 
correlation coefficient for the solid line is R 2 — 0.9987. 



borhood of unity which is important. Indeed, in both 
Figs. [3] and [5] the local linearity in the neighborhood of 
the physical curvature potential is implied. 

The domination of secondary waves for the ERN tails 
for initial data that include a static moment is illustrated 
in Fig. [6l that shows the field as a function of A for a 
number of evaluation times at = 0. Here, the toy 
potential is as in (|18[) . except that the first equation in 
(fT5|) is taken for rv(ro) > > — r*(ro) and the second 
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for |r*| > r*(ro). As shown in Fig. [6j at early evaluation 
times the field is non-linear in A, and at later times it 
becomes linear in A. Also, for A ~ 1, the field is effec- 
tively linear already at earlier times. Also for the case 
of an initially static moment present in the initial data, 
the conclusion is that the tails in ERN are governed by 
secondary waves. Again, also for early times (in the tail 
regime) the local linearity with A in the neighborhood of 
the physical curvature potential is implied by Fig. [6] 



-7.5 - 




0.5 1 1.5 2 2.5 3 3.5 4 

A 

FIG. 6: The magnitude of the tail in ERN at a fixed eval- 
uation point as a function of A for different values of t, for 
the same choice of ro. Same initial data and evaluation point 
(r* = 0) were used in all cases. Here, I = 1, and ro = 100M. 
The evaluation times are: *: t = 300M; +: t = 400M; □: 
t = 500M; o: t = 600M; X: t = 800M; and A: t = 1000M. 

An important result is that for all A ^ 0, the tail 
power-law index remains unchanged, ft is only the am- 
plitude of the tail that is sensitive to A, not the decay 
rate. Same conclusions are found also for the case of 
initial data that include an initially static moment. 



IV. ERN WITH INITIAL DATA OF COMPACT 
SUPPORT 

In this Section we adapt the GPP argument to the case 
of ERN, when the initial data have compact support. 
In Section [V] we consider the case of initial data which 
include an initially static moment. 

Motivated by the preceding considerations, we define 
spacetime duality as follows: 

Definition: Two spacetimes (fixed background geome- 
try and an intial data set for linearized perturbations) 
are said to be dual if the asymptotic contributions of sec- 
ondary waves to the late-time tails in the two spacetimes 
(at a fixed observation point) have the same decay rate. 

Consider an ERN spacetime, with initial data of com- 
pact support. The tails observed on r = const at time 
t have contributions "coming from the right" (incoming 



waves) and contributions "coming from the left" (outgo- 
ing waves) (see the diagram on the left of the top row 
of Fig. [7]). (We first consider the observation point to 
coincide with the maximum of the effective potential for 
reasons to be discussed below, but then our results for 
the tail are independent of the evaluation point such that 
this assumption does not jeopardize the generality of our 
discussion.) Because these are linear waves, we may sep- 
arate the waves on r = const at time t into the two types 
of waves, incoming and outgoing, which make them. We 
next argue, that the tails of the original problem (i.e., 
ERN with perturbations of compact support, and waves 
which are a combination of outgoing and incoming waves) 
are the same as the tails of two superposed problems: 
ERN with the same initial data but only waves coming 
from the right, and ERN with the same initial data but 
only waves coming from the left (right hand side of the 
top row in Fig. [7|). Such a decomposition can be done 
because of two reasons: (i) the linearity of the problem, 
and (ii) our previous result that it is only the secondary 
waves which are important for the generation of the tails. 
(It is important that we discuss only secondary waves. If 
higher-order waves were allowed, then there would be 
no clear separation of waves coming from the right and 
waves coming from the left, as the two would be coupled 
in an intricate manner.) 

Next, consider the case of ERN with the waves coming 
from the right. We argue that the (partial) tails in this 
case are the same as the tails in Schwarzschild for the 
same initial data, because at great distances the effective 
potential in Schwarzschild is the same as in ERN. [De- 
viations are only at 0(r~ 3 ).] Therefore, we can replace 
the first diagram on the right hand side of the top row of 
Fig. [7] with the first diagram on the second row. Similarly, 
the (partial) tails coming from the left in ERN (second 
diagram on right hand side of the top row) are the same 
as the tails in Schwarzschild, if we reflect the initial data 
to the other side of the potential barrier, and change the 
initial data to be initially incoming (second diagram on 
the second row). This is the case because the effective 
potential in ERN close to the event horizon is the same 
as the effective potential in ERN at great distances, and 
the latter in its turn is the same as the effective poten- 
tial at great distances in Schwarzschild. Consequently, 
the ERN effective potential close to the event horizon 
is asymptotically the same as the effective potential at 
great distances in Schwarzschild. It is important that 
we reflect the initial data because we need to keep the 
property of the original diagram, that the direction of 
the waves is the same as the original direction of the ini- 
tial pulse. (The wave propagation is strictly speaking not 
the same in these three spacetimes, because the effective 
potentials at finite distances are not identical. However, 
the tails are generated only by the asymptotic parts of 
the curvature potential, and these are the same in these 
three spacetimes.) 

The original tail problem in ERN is then dual to the 
suporposition of two Schwarzschild problems, which, be- 
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cause of the linearity of the problem, we can recombine 
into a single problem by adding together the original ini- 
tial data and the reflected initial data. There is one more 
subtle point, though: The tails we obtained are not the 
full tails, because in the recombined diagrams we only 
have waves coming from the right, whereas in the full pic- 
ture in Schwarzschild we would have waves coming both 
from the right and from the left. Adding the waves in 
Schwarschild which come from the left, we argue, would 
not change the tails. Allowing for waves coming from the 
left in Schwarzschild would add only negligible amount of 
secondary waves, because in Schwarzschild the effective 
potential drops off exponentially near the event horizon. 
In addition, because the original pulse has compact sup- 
port, there will be no primary waves coming from the 
left at late times. Therefore, we argue that we can add 
to our diagram also the waves coming form the left in 
Schwarzschild, because they will not change the tails. 
This way, we obtain the last diagram in Fig. [7] (third 
row), which is the full diagram for a tail problem in 
Schwarzschild for initial data different from the original 
one (the original initial data superposed with their reflec- 
tion) . As the tail problem in Schwarzschild is completely 
understood (the exponent of the tails is 2£ + 3), we argue 
that the same exponent for the tails is expected also in 
ERN. Below, in Section [AT] we show, using numerical 
simulations, that this is indeed the case. 





FIG. 7: The duality transformations of the tail problem. See 
the text for details. 



V. ERN WITH INITIAL DATA WITH INITIAL 
STATIC MOMENT 

The arguments in this case follow closely those of the 
preceding Section, and will therefore be described here 
briefly. For initial data with an "initially static moment" 
(i.e., a multipolc moment of the static solution is present 
on the initial data hypersurface) , for Schwarschild 



if, 



Sch 



Qi 



M 



M 



1 



J+ 1 



(r > M) 



where Qi is the Legendre function of the second kind, 
while for ERN @] 



^ERN 



(r - M/+ 1 



and 



if, 



ERN 



= (r - My 



(r > M) 



(0 < r - M < M) . 



Recall that the form of the initial data with an initially 
static moment is important only for |r»| 3> M. That is, 
it is only the asymptotic drop off rate of the field (on the 
initial hypersurface) away from the peak of the effective 
potential that determines whether an initially static mo- 
ment is present or not. Remarkably, the initially static 
moment initial data are the same for ERN as they are 
for Schwarzschild for r* ~S> M . We may therefore take 
the initial data in ERN to be V+ RN for r* > and ip ERN 
for r* < 0. Notice that t/j^RN is consistent with the dy- 
namical requirement that a local observer sees a regular 
field as she is crossingthe event horizon. This require- 
ment is presented in [7j as tp ~ const + const/it for all 
£, which at t = becomes ip ~ const + const/|r*|. This 
requirement yields, in fact, the slowest drop off of the 
field. Faster drop off is not disallowed, and will lead to 
a vanishing field on the EH, instead of a non-vanishing 
constant. Notice, however, that the dynamical require- 
ment is consistent with our initial data (at t = 0) only 
for £ = 0, 1, such that for £ > 2 it does not represent our 
choice for initial data. 

In Figure [5] we describe the duality transformations as 
follows. First, the initial value problem, being a linear 
one, is equivalent to the superposition of two problems: 
in the first, ip = V>+ RN for r* > and if> = for r* < 0, 
and in the second, ip — for r* > and ip = ip EKN 
for < 0. At the next step, because of the asymptotic 
symmetry of the ERN effective potential, the tail problem 
in ERN in the second problem above (if) = for r* > 
and ip = t/^RN f or r ^ < q-j j g dual to an ERN problem 

with ip = -0E RN for r* > and ip = for r* < 0. Because 
of the linearity of the scattering problem, we may now 
combine the two problems into a single one, for which 
ijj = ip EKN + ipf BN for r* > and ip = for r* < 0. As for 
r* > M ip E ™ > ipf RN 1 the initial data are dominated 
at large distances by i/^ RN , so that the tail problem is 
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dual to ERN with ip = ^_ for r* > and -0 = for 
r* < 0. Lastly, because of the duality of the tail problem 
in ERN and in Schwarzschild, the original tail problem 
is dual to a Schwarzschild tail problem, where the initial 
data are given by ip ~ r~ e for r* > and ip = for 
r* < 0. 



ERN ERN ERN 




FIG. 8: The duality transformations of the tail problem. See 
the text for details. 

The last tail problem is not a trivial one. It has initial 
data that are stronger than that of an initially static 
moment, but it is similar to the initially static moment 
in having non-compact initial data. A numerical solution 
of this problem, that is Schwarzschild with ip ~ r~ l for 
r* > and ip — for r* < is presented in Fig. [^ which 
displays the local power index n(t) := —tip/ip Q as a 
function of time for three values of I. 
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FIG. 9: The local power index n(t) as a function of m/t in 
Schwarzschild, for initial data given by ip ~ r* for r* > 
and ip — for r* < 0, for a scalar field (s = 0), for I = 1 
(upper panel), £ = 2 (middle panel), and I = 3 (lower panel). 
See the text for details. 



VI. NUMERICAL SIMULATIONS OF TAILS IN 

ERN 

The preceding discussion provides us with the expec- 
tation that initial data in ERN with an initial static mo- 
ment present lead to the same indices of power-law tails 
as in Schwarzschild, namely at late times ip ~ t~( 2e+2 \ 
To test whether these expectations are indeed realized, 
we present in Fig. [TO] the field along r = const as a func- 
tion of time for initial data of an initial static moment I, 
for various values of the latter. In Fig. [TT]we present the 
local power index n(t) for the same data. The numerical 
results agree with our prediction: the late time tails in 
ERN for initial data that are those of an initial static mo- 
ment drop off with time as tp ~ £~( 2 ^+ 2 ) . Finally, in Table 
Uwe confront our prediction for the tail power-law index 
in ERN with numerical results. The numerical results for 
the power-law indices are extrapolations to t oo of the 
local power indices n(t). 
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FIG. 10: The field along r* = as a function of t, for initial 
data with an initially static moment, for a scalar field (s = 0), 
for L = 1 (upper panel) , £ = 3 (middle panel) , and I = 4 (lower 
panel). See the text for details. 

The effective equivalence of the event horizon and fu- 
ture null infinity is further demonstrated by the tails 
along outgoing and incoming null rays at large values 
of advanced and retarded times, respectively. While the 
tails along a r* = const curve fall off at late times 
as H^+H-i) f or Schwarzschild and ERN, where 

/i = 1,2 depending whether an initially static mo- 
ment is present or not, along future null infinity the 
Schwarzschild tails fall off as u~^ l+ ^\ where u = t — r* 
is retarded time. We find this result to remain un- 
changed also for ERN (Fig. [TOJ. The main difference 
between Schwarzschild tails and ERN tails occurs along 
the event horizon: while in Schwarzschild the tails fall off 
as w-^+M-i^ i n ERN they fall off as v~^ +fl \ the same 
as along future null infinity (Fig. [TO)) . Here, v = t + r* is 
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FIG. 11: The local power index n along r, = as a function 
of M/t, for initial data with an initially static moment, for a 
scalar field (s = 0), for I = 1 (left upper panel), £ — 2 (right 
upper panel), £ — 3 (left lower panel), and £ = 4 (right lower 
panel). See the text for details. 



TABLE I: Confrontation of the prediction of this paper based 
on duality arguments (21 + 2) for the power-law index along 
r* = const with our numerical results for ERN and initial data 
of an initially static moment. The relative error is computed 
as the difference between the numerical result and 2£ + 2, 
divided by the latter. 
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FIG. 12: The field along r* = (solid curve), along the event 
horizon (dotted curve), and along future null infinity (dashed 
curve). Upper panel: Initial data of a pure initially static 
moment. Lower panel: Initial data of an initially incoming 
pulse of compact support between — 10M < r, < 10M. For 
both cases 1 = 1. The "bending" of the curves along the 
event horizon and along future null infinity for late advanced 
and retarded times, correspondingly, occurs because these two 
asymptotic regions of spacetimes are represented by large (but 
finite) retarded and advanced times, respectively. 



APPENDIX A: NUMERICAL CODE AND 
CONVERGENCE TESTS 

We have used two versions of the code, that allow us 
to use either Cauchy or characteristic data. For the lat- 
ter, we used a standard characteristic code in 1+1D in 
double-null coordinates u,v. To use Cauchy data, we 
modified the code so that initial data (for ip and for ip) 
are specified at t = 0, but each computational cell is cal- 
culated using the regular characteristic diamond. The 
computational domain is the domain of influence of the 
initial data (characteristic or Cauchy), so that no bound- 
ary conditions are specified. 

In what follows we denote ips — ^{u, v), ips — ip(u,v + 
Av), tpw = tp(u + Au,v), and ipN = tp(u + Au,v + 
Av) . In each computational cell we use the second-order 
algorithm 



IpN = IpE + IpW 



ips ~ j V ipo Au Av . 



In practice, we take = (ipE + ipw + ipN+ips)/^, so that 
we still need to solve for ipN- Collecting the coefficients 
of the fours field points, one finds that 



1p N = 



1 - ^ V Au Av 
1 + ^ V a Au Av 



(ip E +ipw)-ips- (Al) 



This scheme is second order. 

Characteristic data can now be evolved by straight- 
forward marching. Cauchy data presents us with the 
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problem of evolving the initial time step, as no field are 
specified on ipg. To determine the first time step accord- 
ing to the given initial data, consider first the case of 
momentarily stationary initial data. In that case, "0 = 
at t = 0, or ippf = ips on the initial slice. Substituting 
into (|A1|) . we find for momentarily stationary initial data 

1 1 - 4 Vo Au Av , s , , 

V^V = ; , \ 6 T " - — We + ifor). (A2) 



2 1 + i v AuAv 



Two other special cases to consider are tp lV = and ip jU = 
0. In the former case, ipE = ips on the initial slice, which 
we write as ipg = —ipw + ipN + 4>e, so that 



ipw - te v o AuAvipt 



1 + ^ Vo Am Av 



(A3) 



Similarly, if ip iU — 0, one may take ipw = ips on the 
initial slice, or ips = ipw + ipN — ipE, so that 



1p N = 



ipE - jqVqAu Av ip w 



1 + i V Au Av 



(A4) 



General initial data for ip can now be obtained by a linear 



combination of (|A3[) and (|A4[) . In particular, the momen- 
tarily stationary initial data are obtained from (|A3[) and 
(|A4[) by averaging them. 

We have done numerous convergence tests. First, we 
checked our code for a known exact solution (static so- 
lution in Schwarzschild). Then, we tested the global 
convergence order by finding the behavior of vector L p 
norms, and by monitoring the convergence order globally, 
throughout the entire computational domain. As an il- 
lustration for the convergence tests we have done, Fig. [13] 
shows the local convergence order along a = const 
worldline as a function of time. For grid spacings of 
h, ah, a 2 h we calculate the convergence order as 



N = \og a 



Ipah ~ Iph 



Ipah ~ 1pa?h 



where each of the three fields is evaluated at the same 
physical grid point. In all cases we found global second- 
order convergence throughout the computational do- 
main. 
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FIG. 13: Convergence test for a scalar field (s = 0) and dipole 
mode (£ — 1), for Schwarzschild (solid) and for ERN (dashed) 
for initial data with an initial static moment. Shown are 
the convergence order N as a function of time along r» = 
0. Oscillations of the convergence order at the quasi-normal 
mode epoch are typical of oscillatory data. 



